An ab initio study of 3s core-level x-ray photoemission spectra in transition metals 
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We calculate the 3s- and 4s-core-level x-ray photoemission spectroscopy (XPS) spectra in the 
ferromagnetic and nonmagnetic transition metals by developing an ab initio method. We obtain 
the spectra exhibiting the characteristic shapes as a function of binding energy in good agreement 
with experimental observations. The spectral shapes are strikingly different between the majority 
spin channel and the minority spin channel for ferromagnetic metals Ni, Co, and Fe, that is, large 
intensities appear in the higher binding energy side of the main peak (satellite) in the majority spin 
channel. Such satellite or shoulder intensities are also obtained for nonmagnetic metals V and Ru. 
These behaviors are elucidated in terms of the change of the one-electron states induced by the 
core-hole potential. 
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PACS numbers: 79.60.-i 71.15.Qe 71.20.Be 



I. INTRODUCTION 

X-ray spectroscopy has been extensively used for 
studying electronic properties in solids. Core-level spec- 
troscopy is particularly useful for investigating the elec- 
tronic states through the dynamical response to the 
photo-created core-hole. It is well known that the re- 
sponse function in metallic systems exhibits the singu- 
lar behavior near the Fermi edgeii"— In the case of the 
core- level x-ray photoemission spectroscopy (XPS), the 
spectra display the asymmetric shapes as a function of 
binding energy in the vicinity of the threshold 4 Apart 
from the edge singularity, some structures have been ob- 
served in the high binding-energy region in some ferro- 
magnetic transition metals and their compounds. A no- 
table example is a satellite peak on the 2p XPS in fer- 
romagnetic metal Ni, which is located around the region 
6 eV higher from the thresholdi^ Feldkamp and Davis^ 
analyzed these XPS spectra by evaluating the overlap 
between the excited states and the ground state, using 
a numerical method on the linear-combination-atomic- 
orbital model. They clarified the origin of satellite as a 
combined effect of the core-hole screening and the inter- 
action between electrons. 

As regards the 3s XPS, many experiments have al- 
ready been carried out on ferromagnetic transition metals 
and their compounds. ^^^^ Several materials show satel- 
lite or shoulder structures as a function of binding energy, 
which are interpreted as a result of the 3s level splitting 
due to the exchange interaction between the 3s electrons 
and the valence electrons in the polarized 3(i states. In 
some cases, they have been related to the local magnetic 
moment On the other hand, having intensively inves- 
tigated Fe 3s XPS in various iron compounds. Acker et 
al.^ revealed that only poor correlation exists between 
the satellite structures and the magnetic moments. They 
also found that the Fe 3s XPS spectra show the satellite 
structure even in some Pauli paramagnetic compounds. 
Furthermore, having investigated the Mn and Fe 3s XPS 
spectra in the insulating compounds. Oh et alii^ con- 
cluded that the splitting between and main and satellite 



peaks dose not reflect the 2>d moment when the effect of 
the charge-transfer becomes important. 

We have developed an ah initio method to calculate 
the XPS spectra by extending the theory of Feldkamp 
and DaviSi^-.ii2, Here we briefly summarize the procedure 
of calculation. First, we carry out the band structure 
calculation within the local density functional approxi- 
mation (LDA) to obtain the one-electron states in the 
ground state. Next, instead of considering the system 
with only one core hole in crystal, we consider a sys- 
tem of super-cells with one core-hole per cell. Strictly 
speaking, we should consider the former system for the 
XPS event, but the latter system is expected to work bet- 
ter as increasing the cell size. These systems correspond 
to a kind of impurity problem, where the local charge 
neutrality has to be satisfied according to the Friedel 
sum rulei^S We carry out the band structure calculation 
based on the LDA, in which the exchange and Coulomb 
interactions between the core electrons and the conduc- 
tion electrons and between the conduction electrons are 
taken into account through the exchange-correlation po- 
tential. The charge variation due to screening the core 
hole in the final states is also taken into account within 
the super-cell approximation. To guarantee the charge 
neutrality, we add one extra conduction electron in each 
super-cell, and seek the self-consistent solution. Then, 
with the calculated one-electron states, we discretize the 
momentum space into finite number of points, and con- 
struct final states by distributing electrons on these one- 
electron states. The final state with the lowest excited 
energy is given by piling the same number of electrons 
into low energy one-electron levels at each k-point as that 
in the ground state. We prepare the other final states by 
creating one electron-hole (e-h) pair, two e-h pairs, and 
so on. Finally, we calculate the XPS spectra by evaluat- 
ing the overlaps between thus obtained final states and 
the ground state with the help of the one-electron wave 
functions. 

The purpose of this paper is to systematically clarify 
the relation between the spectral shapes and the screen- 
ing process by calculating the spin resolved 3s XPS spec- 



2 



tra in a series of ferromagnetic metals Ni, Co, and Fe. 
The usefulness of our ab initio method is demonstrated. 
We have already reported the spectra in ferromagnetic Fe 
in Ref. [l^. In these metals, spectral shapes have char- 
acteristic dependence on elements and spin channels; the 
spectra have satellite or shoulder in the majority spin 
channel, while the spectra show single peak structures in 
the minority spin channel. Here we define the majority 
(minority) spin channel by the process that the 3s-core 
electron is photo-excited to the vacuum state with the 
same spin as the majority (minority) spin in the con- 
duction band states. As far as we know, such spectra 
have been analyzed only by using a single band Hub- 
bard model, and has been related to the 3s-3d exchange 
interaction.^^ However, the model is, we think, too sim- 
ple to compare the calculated results quantitatively with 
the experimental data and to draw definite conclusion. 

Applying the ab initio method, we calculate the spec- 
tra in good agreement with the experimental observa- 
tions. The screening effects are quite different between 
the spin channels due to the exchange interaction be- 
tween the 3d electrons and the core hole. The 3d states 
are modified by the core-hole potential at the core-hole 
site, and sometimes quasi-bound states are created near 
the bottom of the 3d band. The e-h pair excitations 
from such quasi-bound states to the empty states cor- 
respond to the satellite or shoulder intensities. We find 
that the presence of the quasi-bound states is not sufR- 
cient and the 3c? bands have to be partially occupied in 
the ground state, in order that the satellite or shoulder 
structure appears. These considerations well explain the 
characteristics of the XPS spectra. 

Furthermore, to clearly show that the presence of satel- 
lite or shoulder has no direct relation to magnetic states, 
we calculate the 3s XPS spectra in nonmagnetic metals 
V and Cu, and the 4s XPS spectra in nonmagnetic metal 
Ru. We obtain shoulders in the high binding-energy re- 
gion in V and Ru. On the other hand, we have no such 
structure in Cu, although the localized bound states are 
clearly created below the bottom of the 3d band. We 
could explain these behaviors in the same way as in the 
ferromagnetic metals. 

The present paper is organized as follows. In Sec. II, 
we formulate the XPS spectra with the ab initio method. 
In Sec. Ill, we present the calculated XPS spectra and 
discuss the behavior. The last section is devoted to the 
concluding remarks. 



II. PROCEDURE OF CALCULATION 

A. Formula for XPS spectra 

We consider the situation that a core electron is excited 
to a high energy state with energy e by absorbing an x-ray 
photon with energy ujq and that the interaction between 
the escaping photo-electron and the other electrons could 
be neglected. The probability of finding a photo-electron 



with energy e and spin a could be proportional to 



27T\w\^Y.\(f\'-\9)\'Siuj, + Eg-e- Ej), (I) 
/ 

where w represents the transition matrix element from 
the core state localized at a particular site to the state 
of photo-electron, and is assumed to be independent of 
energy e and spin a. The Sa is the annihilation opera- 
tor of a relevant core electron, which is assumed to have 
only spin a as the internal degrees of freedom. The kets 
\g) and |/) represent the ground state with energy Eg 
and the final state with energy Ef, respectively. We de- 
fine I/) by excluding the photo-electron. In the following 
calculation, we replace the (5-function by the Lorentzian 
function with the full width of half maximum (FWHM) 
2Vs with Ts = 1.0 eV in order to take account of the 
life-time broadening of the core level. 

B. Construction of final and initial states 

In order to simulate the photo-excited states, we con- 
sider a periodic array of super-cells with one core-hole 
per cell, and calculate the one-electron states by means 
of the band structure calculation based on the full po- 
tential linear augmented plane wave (FLAPW) method. 
We use the 3x3x3 bcc super-cell for Fe and V as shown 
in Fig. 1 in Ref. [l^, and the 3 x 3 x 3 fee super-cell 
for Co, Ni, Cu and Ru as shown in Fig. [1] where the 
core-hole sites form a bcc lattice and an fee lattice, re- 
spectively. The larger the unit cell size is, the better 
results are expected to come out. The 3s- or 4s-core 
states in transition metals are treated as localized states 
within a muffin-tin sphere, so that we could specify the 
core-hole site. To ensure the charge neutrality, we as- 
sume He + 1 band electrons per unit cell instead of 
band electrons, where is the number of band electrons 
per cell in the ground state. One additional electron per 
unit cell would not cause large errors in evaluating one- 
electron states in the limit of large unit-cell size. The 
self-consistent potential is obtained as the potential for 
the fully relaxed (screened) state. We write the resulting 
one-electron state with energy eigenvalue ecm(k) as 

V'<Tnk(r) = ^ 0<Tnk(r - R-i) exp(ik • R^), (2) 

with 0crnk(i") = -"crrik (r) e'*' '", where Uo-nk(r) has the pe- 
riod of the super-cell, and j runs over Nc super-cells. 
Needless to say, wave vector k's have Nc discrete val- 
ues in the irreducible Brillouin zone. We use these one- 
electron states as substitutes of the states under a single 
core-hole. We distribute rig band electrons per super- 
cell on these states to construct the excited states. In 
addition, we carry out the band calculation in the ab- 
sence of the core-hole with assuming rig band electrons 
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per super-cell. The wave function and energy eigenvalue 

are denoted as V'i'ik('") ^^'^ etn(k), respectively. All the 
lowest Ne — Nc X He Icvcls are occupied in the ground 
state. 

The final states |/)'s are constructed by using the one- 
electron states calculated in the presence of the core-hole 
in accordance with the following procedure. Defining 
n^^ (k) by the number of levels with spin a and wave 
vector k below the Fermi level in the ground state, we 
distribute ni°^ (k) electrons with spin a and wave vec- 
tor k in the states given in the presence of core hole. 
The final state | /o) containing no e-h pair is constructed 
by distributing electrons from the lowest energy level up 
to the ni°'' (k)'th level with spin a for each wave vec- 
tor k. The final states |/i/)'s containing i/ e-h pairs are 
constructed by annihilating ly electrons in the occupied 
conduction states and creating v electrons in the unoc- 
cupied conduction states from |/o). 




# Core Hole Site 



Figure 1: Sketch of a super-cell containing core- holes in fee 
structure. Core-hole sites indicated by solid circles are as- 
sumed forming a 3 x 3 x 3 fee lattice. The super-cell for the 
bee structure is shown in fig. 1 in Ref. [Ygj . 



C. Overlap integrals 

We assume that the transition matrix elements be- 
tween the core state and the photo-excited states are 
constant. The remaining matrix elements connecting the 
ground and final states are expressed by 



with 
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where the integral is carried out within a unit cell. Sub- 
scripts i = (cr, n, k) and i' = (cr',n',k') are running over 



occupied conduction states in the presence of core hole 
and in the ground state, respectively. We eliminate the 
overlaps between the wave functions for the core levels. 
The corresponding energies difference is given by 

AE = Ef^ - Eg 

^Ef,^-Ef,+Ef,~Eg, (5) 

where Efg — Eg includes the energy of core hole, and is 
treated as an adjustable parameter in the present study 
such that the threshold of XPS spectra coincides with 
the experimental value. The excitation energy Ef^ — Ef^ 
with u e-h pairs is given by 

Eu-Efo = J2i'j-e.), (6) 

(ij) 

where e^'s are the Kohn-Sham eigenvalues, and — e; 
stands for the energy of e-h pair of an electron at level j 
and a hole at level i. Although the Kolm-Sham eigenval- 
ues may not be proper quasi-particle energies, they prac- 
tically give a good approximation to quasi-particle ener- 
gies, except for the fundamental energy gap.'^^-^'^ Substi- 
tuting Eqs. ^ and (O into Eq. H]), we obtain the XPS 
spectra. 

In the actual calculation, instead of Nc k-points, we 
pick up only the F point as the sample states for cal- 
culating XPS spectra. For Ni, we pick up the X point 
(and the equivalent Y and Z points) in addition to the 
F point, since the 3c? band states at the F point are fully 
occupied by both up-spin and down-spin electrons even 
though the 3 x 3 x 3 fee super-cell is used. 

Before closing this section, we briefly mention the XPS 
intensity at the energy of threshold. The final state |/o) 
with the lowest energy (no e-h pair) has a finite overlap 
with the ground state \g), giving rise to intensities at the 
threshold. In principle, such overlap converges at zero 
with Ne ^ oo, according to the Anderson orthogonal- 
ity theorem.^ In such infinite systems, energy levels be- 
come continuous near the Fermi level and thereby infinite 
numbers of e-h pairs could be created with infinitesimal 
excitation energies, leading to the so called Fermi edge 
singularity in the XPS spectra. The finite contribution 
obtained above arises from the discreteness of energy lev- 
els and could be interpreted as the integrated intensity 
of singular spectra near the threshold, in consistent with 
the model calculations for other systemsi^i^i 



III. RESULTS AND DISCUSSION 

A. Ferromagnetic Transition Metals 

In this subsections, we refer to majority(minority) spin 
as up(down)-spin. The Ni metal takes an fee structure. 
For simplicity, the Co metal is assumed to take an fee 
structure, although it actually takes an hep structure. 
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Table I: Screening electron number in the d-symmetric states 
inside the muffin-tin sphere at the 3s core-hole site. The radii 
of the muffin-tin spheres are 2.0 Bohr. 

3s hole spin Audf Anai Audf + Aridi 



Fe 


up 


-1.44 


2.38 


0.94 




dn 


0.47 


0.48 


0.95 


Co 


up 


-0.52 


1.55 


1.03 




dn 


0.24 


0.75 


0.99 


Ni 


up 


0.07 


0.87 


0.94 




dn 


0.26 


0.65 


0.91 



The Fe metal takes a bcc structure. Figures [5] and [3] 
show DOS's projected onto the states with d symmetry 
(d-DOS) at the core-hole site for Ni and Co, respectively. 
The corresponding DOS's for Fe are shown in Fig. 3 in 
Ref. [l^. In these calculations, six fc-points are picked 
up in the irreducible Brillouin zone for super-cell sys- 
tems. The DOS's calculated with no core-hole are essen- 
tially the same as those reported by Moruzzi, Janak and 
Williams. Table U lists the screening electron number 
Anda hi the d-symmetric states with spin a, that is, the 
difference of the occupied electron number between in the 
presence and in the absence of the core hole inside the 
muffin-tin sphere. 
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On the basis of these one-electron states, we calcu- 
late the 3s XPS spectra, by following the procedure de- 
scribe in Sec. II. Figures SI [S] and IH] are the spectra thus 
calculated as a function of the binding energy Uq — e 
for Ni, Co, and Fe, respectively, in comparison with the 



experiments J 



The spectral shape in Fig. [6] for Fe is 



§htly different from our previous result (Fig. 4 in Ref. 
, since the present calculation takes full account of ex- 
citations up to three e-h pairs in comparison with only 
up to two e-h pairs in Ref. [Tol The spectra are strikingly 
different between the up-spin channel and the down-spin 
channel and strongly depend on elements, in good agree- 
ment with the experiments. In the following, we explain 
the origin of these behaviors in relation to one-electron 
states screening the core hole. 



Figure 2: Calculated d-DOS at the core-hole site in the super- 
cell system in ferromagnetic nickel; (a) d-DOS with no core- 
hole, (b) d-DOS when the 3s up-spin electron is removed, (c) 
d-DOS when the 3s down-spin electron is removed. 
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Figure 3: Calculated d-DOS at the core-hole site in the super- 
cell system in ferromagnetic cobalt; (a) d-DOS with no core- 
hole, (b) d-DOS when the 3s up-spin electron is removed, (c) 
d-DOS when the 3s down-spin electron is removed. 



1. Up-spin channel 

First we consider the situation that a up-spin 3s elec- 
tron is removed in a unit cell. As shown in Figs. [2] (b) 
and[3l^b), the d-DOS at the core- hole site are strongly 
modified by the core-hole potential for Ni and Co. The 
situation is similar to Fe, as shown in Fig. 3 in Ref. [l^. 
The d-DOS's of the down-spin states are stroirgly pulled 
down, forming quasi-bound states around the bottom of 



Table II: Absolute squares \Af\^ and |j4^|^, where Af and A^ 
represent the up- and down-spin parts of the overlap integral, 
respectively, between |/o) and Scr\g} in the cr-spin channel, 
that is, {fo\s:y\g) = A^Ai. 





up-spin core-hole 


dn-spin core-hole 








\A,f 




Ni 


0.940 


0.234 


0.770 


0.572 


Co 


0.983 


0.049 


0.980 


0.601 


Fe 


0.967 


0.127 


0.969 


0.909 



the 3d band. On the other hand, the d-DOS's of the up- 
spin states are slightly pulled down for Ni, and pushed 
upward to the higher energy region for Co and Fe. The 
d-DOS's at the site without a core-hole are essentially 
the same as those in the ground state with no core-hole 
Since up-spin electrons are prevented from coming 
close to the core-hole site due to the exchange interac- 
tion with the up-spin core hole, the screening is almost 
done by down-spin electrons. This tendency is clear in 
Ni; Audi = 0.87, while An^t = 0.07, as shown in TableU] 
For Co and Fe, the core-hole potential is overscreened by 
down-spin electrons; Audi — 1.55(Co), Aridi — 2.38(Fe). 
This overscreeniirg is compensated by up-spin electrons; 
Aud^ = -0.52(Co), Aud^ = -1.44(Fe). As a resuh, the 
screening electron numbers become almost unity: 0.94, 
1.03, and 0.94 for Ni, Co, and Fe, respectively, indi- 
cating that the screening is nearly completed by the d- 
symmetric states. The magnetic moments at the core- 
hole site are — 0.2/iB, — 0.5/iB, and — 1.8/j,b for Ni, Co, 
and Fe, respectively, which are opposite to those without 
tlie core-hole, 0.6 /ie, 1.6/j,B) 2.1 /ib, for Ni, Co, and Fe, 
respectively. 

We note that the magnitude of the overscreeniirg by the 
down-spin electrons becomes weaker Fe, Co, and Ni in 
that order. In Ni, the number of the down-spin electrons 
screening the core-hole is just about unity. Consequently, 
the up-spin 3d states become to be less pushed upward to 
the high energy region in that order. In Ni, the effects of 
the core-hole potential and the Coulomb repulsion from 
the screening down electrons subtly cancel out each other 
and the d-DOS of up-spin states is hardly modified. 

It is inferred from these changes in the d-DOS's that 
one-electron wave functions are largely modified by the 
core-hole potential particularly for down-spin states. It is 
necessary to use both the occupied and unoccupied states 
of the ground state in order to expand those modified 
one-electron wave functions for the down-spin electrons, 
since the down-spin 3d bands are partially occupied in 
the ground state. For this reason, the absolute square 
1^4.1 becomes rather small, as shown in table|lll Here 
(A^) represents the the up(down)-spin part of the overlap 
integral between the lowest-energy final state |/o) (no e-h 
pair) and s^lff), and thereby {fo\s^\g) = A^A^. 

The one-electron wave functions for up-spin electrons 
are also modified from those in the ground state. In spite 
of such modification, 's are nearly unity as shown in 
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Table HIl This could be understood as follows. Since the 
up-spin 3d bands are almost fully occupied in the ground 
state, up-spin one-electron states constituting the final 
state l/o) could be represented by a unitary transform of 
those constituting the ground state \g). Therefore, since 
the determinant is invariant under unitary transforma- 
tion, A^'s are close to unity. 

Final states |/;y)'s containing up-spin e-h pairs could 
give rise to only small intensities, since the states of the 
excited electrons with up-spin are almost orthogonal to 
occupied states with up-spin in the ground state \g), and 
thereby the overlap determinants would vanish. On the 
other hand, the final states |/i)'s containing one down- 
spin e-h pair could give rise to considerable intensities, 
since the corresponding one-electron wave functions con- 
tain the amplitudes of the unoccupied one-electron states 
in the ground state \g), and thereby the overlap deter- 
minants would not vanish. Considering various combina- 
tions of one e-h pair, we obtain intensities distributed in 
a wide range of binding energy. 
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to the unoccupied states with down-spin. This excitation 
may be considered as a core- hole plus (P, since one of 
the quasi bound states, which are almost localized at the 
core-hole site, is empty. In this calculation, however, the 
states are not split off from the band bottom edge, indi- 
cating that the core-hole-dn-hole states are only weakly 
bound. The satellite binding energy of the calculated 
spectra is 1 eV smaller than that of the observed spec- 
tra. This discrepancy might be owing to the LDA. The 
Excitations of two e-h pairs would give rise to consider- 
able intensities in the energy range ujq — e = 110 '--^ 120 
eV. Excitations of three e-h pairs give rise to only small 
intensities. 
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Figure 5: 3s XPS spectra in ferromagnetic cobalt as a function 
of binding energy, (a) and (b) are for the up-spin and down- 
spin channels, respectively. The experimental data are taken 
from Ref. [l3|. 



Figure 4: 3s XPS spectra in ferromagnetic nickel as a function 
of binding energy, (a) and (b) are for the up-spin and down- 
spin channels, respectively. The experimental data are taken 
from Ref. 

Figure[Ua) shows the calculated spectra for Ni in com- 
parison with experimental observations in the up-spin 
channel. The calculated spectra have maximum intensity 
at the threshold around ojq — e = 110 eV and the signifi- 
cant satellite intensity around ujq — e = 115 eV. Note that 
that the final states containing one down-spin e-h pair 
give rise to considerable intensities, extending over main 
and satellite regions. The satellite intensity corresponds 
to excitations of one e-h pair from the quasi-bound states 



Figure \^a) shows the calculated spectra for Co. The 
calculated spectra have a broad peak structure with the 
maximum intensity at the threshold around ujq~ e — 101 
eV. They also have large shoulder intensities around 
ujq — e = 104 eV. The former peak originates from the 
excitations with one and two e-h pairs. The contribution 
of the lowest-energy final state |/o) (no e-h pair) is quite 
small due to small \Ai\ . The latter shoulder originates 
from excitations of two and three e-h pairs, probably in- 
cluding the excitations from the quasi-bound states to 
the unoccupied states with down-spin. 
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Figure 6: 3s XPS spectra in ferromagnetic iron as a function 
of binding energy, (a) and (b) are for the up-spin and down- 
spin channels, respectively. The experimental data are taken 
from Ref. 

Figure [S^a) shows the calculated spectra for Fe. The 
calculated spectra consist of a peak around ujq ~ e — 92 
eV and a satellite peak around 95 eV. The satellite peak 
is larger than the peak around the threshold. The final 
states l/i) 's containing one down-spin e-h pair give rise to 
the satellite intensity. The final states |/2)'s containing 
two e-h pairs give rise to a shoulder to the satellite around 
= e = 97 ~ 100 eV, as shown in the figure. The 
excitations of three e-h pairs gives rise to finite but small 
intensities in the wide energy range around 98 eV. 



total screening electron numbers are almost unity, 0.91, 
0.99, and 0.95 for Ni, Co, and Fe, respectively. The local 
magnetic moments at the core-hole site are not signifi- 
cantly changed from those in the ground states. 

The |A-|-p's are not far from unity as listed in Table HIl 
by the same reason as in the up-spin channel. Note that 
jAj^P's are larger than in the up-spin channel, indicating 
that one-electron wave functions for down-spin conduc- 
tion electrons are less modified in the down-spin channel 
than in the up spin-channel. 

Figure IHb) shows the calculated spectra for Ni, in 
comparison with the experiment. The final states |/i)'s 
containing one e-h pair with down-spin give rise to con- 
siderable intensities in a wide energy region 110 ^ 115 
eV, while the final states |/2)'s containing two e-h pairs 
give rise to small intensities in a wide energy region with 
the maximum around Uq — e = 115 eV. We obtain a 
asymmetric shape with a tail in the high-energy region 
in agreement with the experiment. 

Figure O^b) shows the calculated spectra for Co. The 
lowest-energy final states |/o) gives rise to a peak at the 
threshold around ujg — e = 100 eV. Different from Ni, 
final states |/i)'s and |/2)'s give rise to intensities on a 
limited region near the threshold. This suggests that one- 
electron wave functions are modified only for levels in the 
vicinity of the Fermi level. Contributions of final states 
|/3)'s (three e-h pairs) are found negligible. 

Figure mb) shows the calculated spectra for Fe. The 
final state |/o) gives rise to a main peak at the thresh- 
old — e = 91 eV with the largest contribution in the 
three metals. Final states |/i)'s with one down-spin e-h 
pair give rise to a broad peak around 94 eV, which is 
quite small in comparison with the main peak. Although 
the observed spectra show small satellite intensity around 
Ljq ~ e = 97.5 eV, the present calculation gives no inten- 
sity there. It is unclear on the mechanism giving the 
intensity to our knowledge. 



2. Down-spin channel 

When a down-spin 3s electron is removed in a unit 
cell, the screening behavior is quite different from the 
situation where a up-spin 3s electron is removed. As 
shown in Figs. [^Tc) andlHJc), and the bottom panel in 
Fig. 3 in Ref. [l^, the d-DOS's at the core-hole site for 
Ni, Co, and Fe are strongly modified by the core-hole 
potential. Different from the up-spin channel, the effect 
is larger for up-spin conduction states than for down-spin 
conduction states; large weights are transferred to the 
bottom of the conduction band in the up-spin d-DOS's, 
while the weights are slightly shifted downward to the 
lower energy region in the down-spin d-DOS's. 

The screening electron numbers are rather smaller in 
the up-spin state than in the down-spin state, as listed in 
Table m This difference arises from the fact that up-spin 
id states are almost occupied in the ground state. The 



B. Nonmagnetic metals 

In ferromagnetic metals Fe, Co, and Ni, the XPS spec- 
tral shapes remarkably depend on the spin channel, ow- 
ing to the exchange interaction between the core hole and 
the conduction electrons. Although the spectral shape is 
found closely related to the filling of band states, the pres- 
ence of the satellite in the XPS spectra seems to have no 
direct relation to ferromagnetic states. In this subsection, 
taking up typical nonmagnetic metals Cu, V and Ru, we 
clarify this issue. Since the spectra are independent of 
the spin channel in nonmagnetic metals, we consider the 
down-spin channel in the following. We assume the fee 
structure for Cu and Ru, and the bcc structure for V. 

Figure [3 show the d-DOS for Cu. The DOS calcu- 
lated with no core-hole are essentially the same as those 
reported by Moruzzi, Janak and Williams;^^ The d-DOS 
at the core-hole site is shifted downward to the deeper en- 
ergy region. The localize bound-states are clearly created 
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below the 3d bands in both up- and down-spin states. 
As hsted in Table Hill the screening electron numbers are 
given by An^t = 0.31, Audi — 0.24, and the total screen- 
ing electron number is given by An^^-I-And^ — 0.55. This 
indicates that the core-hole potential is not sufficiently 
screened by the 3d electrons at the core-hole site. Here we 
note that the screening mechanism seems to be somewhat 
different from the other metals which have partially oc- 
cupied c?-bands and do not show the bound states. Since 
the 3d band states are almost fully occupied in the ground 
state, the change of the charge density is hardly achieved 
by a unitary transform of the 3d band states. Since the 
bound states are split off from the band bottom edge, the 
radial part of the local atomic wavefunctions constituting 
the bound states at the core-hole site is slightly shrunk 
compared to that in the ground state with no core-hole, 
leading to the small change of the charge density. This 
shrink of radial part of the local atomic wavefunctions 
cannot be described by a unitary transform of the 3d 
band states in the ground state. The states with much 
higher energy in the ground state necessarily constitute 
the bound states to some extent. 
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Figure 7: (a) DOS calculated in a system of super-cell with 
no core-hole in nonmagnetic copper. The solid line represents 
the DOS projected onto the d-symmetry within the mufRn-tin 
sphere, (b) DOS at the site of 3s down-spin core-hole. The 
thin line represents the total DOS divided by the number of 
atoms in a unit cell. 



Figure M shows the d-DOS for V. The d-DOS at the 
core-hole site is shifted downward to the deeper energy 
region for up-spin states, while the d-DOS is shifted up- 
ward for down-spin states. No bound state is formed. 
As listed in Table Hill the screening electron numbers are 
given by An^^ = 1.21, An^i — —0.21, and the total num- 
ber by Audf + Andi = 1.00. The core- hole potential is 
overscreened by up-spin electrons, and the overscreening 
is compensated by down-spin electrons. Note that the 
screening is complete within the 3d electrons at the core- 
hole site. The screening is effective because the 3d bands 
are partially occupied. 
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Figure 8: (a) DOS calculated in a system of super-cell with no 
core-hole in nonmagnetic vanadium. The solid line represents 
the DOS projected onto the d-symmetry within the muffin-tin 
sphere, (b) DOS at the site of 3s down-spin core-hole. 
The solid line represents the DOS projected onto the d- 
symmetry within the muffin-tin sphere. 

Figure M shows the d-DOS for Ru. The d-DOS's at 
the core-hole site are shifted downward to the deeper 
energy region with both spin states, but the change is 
the smallest in the three cases without any bound states. 
This is probably related to the fact that the 4d electrons 
are more itinerant than the 3d electrons. As listed in 
Table lllli the screening electron numbers are An^^ — 
0.77, Andi = 0.24, and Audf + Audi = 1.01, indicating 
that the screening is nearly completed by the 4d electrons 
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at the core-hole site. 
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Figure 9: (a) DOS calculated in a system of super-cell with no 
core-hole in nonmagnetic ruthenium. The solid line represents 
the DOS projected onto the d-symmetry within the muffin-tin 
sphere, (b) DOS at the site of 4s down-spin core-hole. 



Table III: Screening electron number with respect to the d 
symmetry within the muffin-tin sphere at the 3s (4s) down- 
spin core-hole site. The radii of the muffin-tin spheres are 2.0, 
2.1, and 2.3 Bohr for Cu, V, and Ru, respectively. 





3s hole spin 






Andt -1- Audi 


Cu 


dn 


0.31 


0.24 


0.55 
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dn 


1.44 


-0.44 


1.00 


Ru 


dn 


0.77 


0.24 


1.01 



We calculate the up- and down-spin parts of overlap 
integrals between the lowest-energy final state |/o) and 
the ground state \g), which values are listed in Table HVl 



Table IV; Absolute squares \Af\ and |^^| , where Af and A_i 
represent the up- and down-spin parts of the overlap integral 
between |/o) and s^|p), that is, (/o|s4,Isf) = A^A^. 



dn-spin core-hole 








Cu 


0.910 


0.935 


V 


0.571 


0.938 


Ru 


0.681 


0.949 



For Cu, and are close to unity, although 

the one-electron wave functions are strongly modified. 
Since the 3d bands are fully occupied in the ground 
state, one-electron wave functions constituting |/o) are 
nearly expressed by a unitary transform of those con- 
stituting \g), except for the shrink of the radial part of 
the atomic wavef unctions at the core-hole site. There- 
fore, the squares of the overlap determinant \A^\'^ and 
|^4.P are nearly unity. Final states |/i) (one e-h pair), 
1/2) (two e-h pairs), and so on, could have merely very 
small overlaps with si\g), since the one-electron state on 
which the excited electron sits is nearly orthogonal to the 
one-electron states constituting the ground state. Thus 
we have a simple single peak structure coming from |/o) 
without any noticeable intensities on the higher binding 
energy side, as shown in Fig. [TOl This result is consis- 
tent with the experimental observation.^ Note that the 
change of the wave functions is not directly related to the 
spectral shape. The effect due to forming strong bound 
states is not seen. 



Cu 3s XPS 



0.3 



total 
lfo> 



0.2- 



B 0.1 
a 



10 5 -5 

Binding Energy (eV) 



-10 



-15 



Figure 10: 3s XPS spectra in nonmagnetic copper as a func- 
tion of binding energy in the down-spin channel. 



For V, lA^l" is close to unity. This suggests that one- 
electron wave functions for down-spin electrons are little 
modified from those in the ground state. On the other 
hand, \A-f\^ is rather smaller than unity. This suggests 
that up-spin one-electron wave functions constituting |/o) 
include the amplitudes of the unoccupied one-electron 
states in the ground state. In such a situation, the final 
states |/i)'s containing one e-h pair with up-spin could 
have finite overlaps with the ground state. Figure [Tl] 
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shows the calculated spectra. We have a main peak com- 
ing from the final state |/o) at the threshold and the 
noticeable shoulder coming from the final states |/i)'s. 
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Figure 11: 3s XPS spectra in nonmagnetic vanadium as a 
function of binding energy in the down-spin channel. 

For Ru, is again close to unity by the same reason 
as in V. The \A^\^ is smaller than unity, although it is 
larger than that in V. This indicates that one-electron 
wave functions constituting |/o) are less modified by the 
core-hole potential in comparison with V. Final states 
|/i)'s containing one e-h pair with up-spin give rise to a 
shoulder structure with a little smaller intensity than in 
V. 
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Figure 12: 4s XPS spectra in nonmagnetic ruthenium as a 
function of binding energy for the 4s down-spin core-hole. 



for the up-spin channel, while the intensities are con- 
centrated near the threshold with no satellite peak for 
the down-spin channel. The origin of such behavior has 
been explained in relation to the 3d band modified by the 
core-hole potential and the overlap integral between the 
final states and the ground state. Bound or quasi-bound 
states are formed by the core-hole potential, and the e-h 
excitations from such quasi-bound states to the unoccu- 
pied levels would usually give rise to satellite intensities. 
However, the presence of the quasi-bound state is not 
a sufficient condition to the presence of satellite; the ri- 
band should be partially occupied in the ground state, 
and thereby the one-electron wave functions constituting 
the final states include the amplitudes of the unoccupied 
one-electron states in the ground state. If the d-band is 
fully occupied, the satellite intensity would not come out 
even in the presence of the bound state. Note that the 
satellite peak position has no direct relation to the 3s 
level exchange splitting; the LDA calculation gives such 
splittings as 0.7, 1.9, and 2.5 eV for Ni, Co, and Fe, re- 
spectively. 

These results indicate that the presence of satellite is 
not directly related to the ferromagnetic ground state. 
We have clarified this point by calculating the spectra in 
nonmagnetic metals Cu, V, and Ru. For V and Ru, we 
have obtained shoulder structures in the XPS spectra, 
although the structure is rather small for Ru. The origin 
of these behaviors is the same as in the ferromagnetic 
metals. For Cu, only a symmetric peak is found with no 
structure, although the localized bound states are clearly 
formed below the bottom of the conduction band. This is 
because the 3d band is completely occupied in the ground 
state. 

We have calculated the XPS spectra in Ni, Co, Fe, 
and Cu in good agreement with the experiment, while 
we could not find experimental XPS data for V and Ru. 
Acker et al. observed the satellite structures even in some 
Pauli paramagnetic Fe compoundSjS The present results 
would provide an interpretation of their findings. Fi- 
nally, as regards the L-edge spectra, experimental data 
for XPS spectra and the x-ray absorption spectra are ac- 
cumulated, and ab initio approach has been triedi^ The 
extension of the present method to the L-edge spectra is 
left in future study. 
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